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First the second-order perturbations of nonzero-A cosmological models are derived explicitly with 
an arbitrary potential function of spatial coordinates, using the nonlinear version of Lifshitz's method 
in the synchronous gauge. Their expression is the generalization (to the nonzero-A case) of second- 
order perturbations in the Einstein-de Sitter model which were derived previously by the present 
author. Next the second-order temperature anisotropies of Cosmic Microwave Background radia- 
tion are derived using the gauge-invariant formula which was given by Mollerach and Matarrese. 
Moreover the corresponding perturbations in the Poisson gauge are derived using the second-order 
gauge transformations formulated by Bruni et al. In the second-order it is found in spite of gauges 
that tensor (gravitational-wave) perturbations and vector (shear) perturbations without vorticity 
are induced from the first-order scalar perturbations. These results will be useful to analyze the 
nonlinear effect of local inhomogeneities on Cosmic Microwave Background anisotropies. 
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I. INTRODUCTION 

The anisotropies of the Cosmic Microwave Background radiation (CMB) gives us an important information of 
the structure of our universe. So far the relation of CMB anisotropies to inhomogeneities along light paths has 
been analyzed using the linear theory of gravitational instability. Recently the measurements of the anisotropies 
has however become more and more precise, so that we can catch some information related to nonlinear effect 
of inhomogeneities such as in the form of non-gaussianity [1]. 

The general-relativistic second-order nonlinear theory of gravitational instability was studied previously by 
the present author [2, 3], in connection with structure formation, by extending Lifshitz's linear theory [4] to 
the second-order smallness. They were restricted to the vanishing A model (Einstein-de Sitter model) and 
the treatment in the synchronous gauge. Recently the nonlinear gauge transformations and the condition of 
second-order gauge-invariantness have been studied by Bruni et al.[5] and the second-order tranformation from 
the synchronous gauge to the Poisson gauge has been performed by Matarrese et al.[6], and the second-order 
temperature anisotropy of CMB radiation has been derived by Mollerach and Matarrese [7] . The second-order 
theory has, moreover, been extended to useful cases [8, 9, 10, 11] of nonzero-A and non-vanishing pressure, to 
analyze non-gaussianity in the CMB anisotropies. 

In this paper we will derive the second-order perturbations corresponding to the first-order scalar pertur- 
bations, in the nonzero-A and pressureless case using the synchronous gauge, in which the perturbations are 
expressed explicitly with an arbitrary potential function (_F) of spatial coordinates, similarly to our previous 
work [2]. 

Next we represent the second-order CMB anisotropies (82T/T) using the derived metric perturbations in 
the synchronous gauge. The CMB anisotropies, which was derived by Mollerach and Matarrese[7] in arbitrary 
gauges, are gauge- invariant, so that our expressions will be useful in the same way as those in the other gauges. 

Lastly we derive the second-order perturbations in the Poisson gauge from those in the synchronous gauge, 
using the nonlinear gauge transformation exploited by Bruni et al.[5] and Matarrese et al.[6] The latter gauge 
may be familiar and comprehensive, because the line-element in the Poisson gauge is similar to that in the 
Newtonian and Post-Newtonian approximation [12]. 

It is found that, in the second-order, tensor (gravitational-wave) perturbations and vector (shear) perturba- 
tions without vorticity are induced from the first-order scalar perturbations. 

In Appendix A the derivation of basic equations in the extended version of Lifshitz theory is reviewed. 
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II. SECOND-ORDER PERTURBATIONS IN THE SYNCHRONOUS GAUGE 

First our background universe is assumed to be described by isotropic and homogeneous pressureless cosmo- 
logical models which are spatially flat, and their spacetimes are expressed by the line-element 

ds 2 = g^dx^dx" = a 2 {rf)[-dr, 2 + (2.1) 

where the Greek and Latin letters denote 0,1,2,3 and 1,2,3, respectively, contrary to the notation in the 
previous paper [2], and Sij(= Sj = S lj ) are the Kronecker delta. The conformal time t](= x°) is related to the 
cosmic time t by dt — a(r])dri. 

In the comoving coordinates, the velocity vector and the energy- momentum tensor of the pressureless matter 
are expressed as 

u° = l/a, u l =0 (2.2) 

and 

T ° = -p, T° = 0, 7? = 0, (2.3) 

where p is the matter density. From the Einstein equations, the equations for p and the scale factor a are 
obtained : 

pa 2 = 3(a'/a) 2 - Aa 2 , (2.4) 

and 

pa 3 = po, (2.5) 

where a prime denotes d/drj, A is the cosmological constant, and po is an integration constant. 

Next we consider first-order perturbations of the scalar type. The perturbations of metric, matter density 
and velocity are represented by Si 9^v{= h^), Si p, and Si u^. 

When we adopt the synchronous coordinates, the metric perturbations satisfy the condition 

/ioo = and h oi = 0. (2.6) 

The scalar-type solutions of the perturbed Einstein equations are classified into the growing case and the 
decaying case, and the remaining components are expressed in both cases as follows: 



(1) the growing case 

h{ - siF+p{r,)F\i, 



Su" = 0, 8u l = 0, 
i i 

= M^'- 1 )^* (2 - 7) 

where F is an arbitrary function of spatial coordinates x x ,x 2 and x 3 , A = V 2 , h\ = g^hu and P(r]) satisfies 



Its solution is expressed as 



P» + ^-p' -1 = 0. (2.8) 



P= [ V dr)'a- 2 (T]') f d V "a 2 (r)"). (2.9) 
Jo Jo 

The three-dimensional covariant derivatives \i are defined in the space with metric dl 2 = 5ijdx l dx^ and their 
suffices are raised and lowered using Sij, so that their derivatives are equal to partial derivatives, i.e. F^ — Fjj, 
where F, = dF/dx' . 
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(2) the decaying case 



H = P(v)F tij , 

Su° = 0, 8u l = 0, 
i i 



where P(rj) satisfies 



and the solution is 



Sp/p = ^P'AF, (2.10) 
l pa z 



P" + ^Lp> = o (2.11) 



P = [" dr]'a- 2 {r)'). (2.12) 



These first-order density perturbations are consistent with the gauge-invariant variable e m (defined by 
Bardecn[14]), which is described by the equation 

e m " + -ej - \pa 2 e m = (2.13) 

CI £ 

in the pressureless case without anisotropic stresses and entropy perturbations. In fact <5i p/p in Eqs.(2.7) and 
(2.10) satisfy Eq. (2.13) for P in Eqs.(2.9) and (2.12). In both cases, the potential function F is determined by 
the last equation of Eq.(2.7) or Eq.(2.10) which is regarded as the cosmological Poisson equation, if the spatial 
distributions of <5i p/p are given at a specified epoch. 

Now let us derive the second-order perturbations 82 9^v{= &nv),82 P, and 82 u M , where the total perturbations 
are 

8u» = Su^ + SuL 1 , 

Sp/p = Sp/p + Sp/p. (2.14) 

1 2 

Here from the synchronous gauge condition, we have 

£ 00 = and l i = 0. (2.15) 

Then, by solving the equations (A16) and (A17) for lij which are derived from the perturbed Einstein equations 
in Appendix A, we obtain the following expressions of 82 9nv{= t-fiv) and using Eqs. (A17) and (A18) we obtain 
82 p and 82 u^. 

(1) the growing case 

The metric perturbations reduce to 

l\ = P{n)L{ + P 2 (jj)Mi + Q{n)N\i + CI (2.16) 
where = S jl N\ H — N^j and Q(r)) satisfies 

On' ^ 

Q" + —Q' = P-^{P'f (2.17) 

CI Zi 

and the solution is expressed as 

Q = £ dia- 2 ^') £ dr/'a 2 ^') [P( V ") \{P'{v")) 2 ] ■ (2-18) 

In the case A = 0, we have P — (5/2)(P') 2 = because of a oc rj 2 and P = ?7 2 /10, so that Q vanishes. The 
functions L\ and M- are defined by 



V 1 

Li ~ 2 



SFiFj - 2F ■ F tij + ^SyFjFj 
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Mf = ^{l9F i; i^ - 12F tij AF - 3% [f m F m - (AF) 2 ] } 



and N is defined by 



AiV 



28 



(AF) 2 - F M F t 



ki 



The velocity and density perturbations are found to be 



<Su° = 0, Su l = 

2 2 



(2.19) 



(2.20) 



(2.21) 



and 



Sp/p = ^{ l -{l--P l ){2>F,F^ + %FAF)+ l -P[{AFf+F M F M ] 
2 2pa z 1 2 a 2 



1 

+ 4 



'{P'f - ~Cr] [(AF) 2 - F M F M ] - l -^p P '[AF M F M + 3(AF) 2 ]}. 
The last term C\ satisfies the wave equation 

UCl^^{P/afGl+ l -[P-\{P') 
where the operator □ is defined by 

2a' 



l\2 



□0 ee = -a' 2 ^/dr? + -^-d/d V - A) , 



(2.22) 



(2.23) 



(2.24) 



for an arbitrary function <ft by use of the four-dimensional covariant derivative ;, and G\ and G\ are expressed 
as 

G\ ee A(F tij AF - F til F tjl ) + (F^F M - F tik F jt ) M - i%A[(AF) 2 - F M F M ], 

G\ ee FyAF - FaFj, - ^%[(AF) 2 - F fcI F fcI ] - 7N ir (2.25) 

These functions satisfy the traceless and transverse relations 

G\ = 0, Gl = 0, 

G\ = 0, G\j = 0, (2.26) 

so that Cf also satifics 



c/ = o, C^ = 0. 



(2.27) 



This means that Cf represents the second-order gravitational radiation emitted by first-order density pertur- 
bations. The solution of the above inhomogeneous wave equation (Eq.(2.23)) can be represented in an explicit 
form using the retarded Green function for the operator □ [2, 16, 17]. 



(2) the decaying case 

The metric perturbations are 



1 = \p 2 { n ){2F, u F s - AF • F 4 , + ^[(AF) 2 - F M F M ]} + C{, 
where P(rj) in this case is given by Eq.(2.12). The last term C\ is described by the wave equation 

VCt = -l(P/a) 2 Gi, 



(2.28) 



(2.29) 



where G\ satisfies Eq.(2.25). 
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The velocity and density perturbations are 

Su° = 0, 5u l = (2.30) 

2 2 

and 

Sp/p = ^P'{p'[(AF) 2 - F M F M ] - ^P[(AF) 2 + 3F M F M }}. (2.31) 



Here 82 p/p is the matter density perturbations observed by comoving observers and C{ represents the gravita- 
tional radiation emitted by the first-order density perturbations. 

By the way we consider the rotational velocity u a and the corresponding scalar quantity uj defined by 

w « = ^ v ^^ upu ^ and w = ( w « Wq )V2. (2.32) 

In the above perturbations, it is clear that uo a vanishes, because Si = S2 = 0, and so w also vanishes. 

Next, for a later use we here express our solutions (in the growing case) using the notation which was employed 
by Matarrese, Mollerach and Bruni [6, 7] for the gauge-invariant treatment of second-order perturbations. In 
their notation our perturbations are expressed in the following form: 



9ij 

and 



.goo - -a 2 ( V )[l + 2^+^% 

a 2 (v){[l 2^) - ^]8 l0 + xg) + i X g)}, (2.33) 



X& ) =4 )II +4 )T , (r = l,2), (2.34) 



where || and T denote the scalar and tensor perturbations, respectively, and 

1 
3 



Dij = didj - hijA. (2.35) 



The velocity and density perturbations are 

= - 

a 

and 



W 4. v » 4. I,,** 

°0 + W (l) T" 2 ^ (2 ) 



(2.36) 



+ 5^)^+ i^( 2 )p. (2.37) 



P = P(o) ■ - r . 2 
For our perturbations in the synchronous gauge, we have 

V-M = = Z V = = (2.38) 
and the other components are expressed by use of our notation as 

<t>™ = -h/&=- l -{F+ l -PAF), 

4> {2) = -l/2> = -\{PL\ + P' 2 Mt + QAN), 

X™ = =h{ - l -5lh = p(F tij ^AF) or X (1)l1 = PF, 

Y (1)T = 

\xf^ = 4 - = p H + p2M i + QNij - \Sij{PLk + P 2 M£ + QAN), 

iy (2)T = C 1 

V (l) = V (2) = °. 

S 1 p = Sp, \5 {2) P = 5p, (2.39) 

1 £ 2 

where h = h\ and I = It is interesting that xjj^" includes not only scalar perturbations, but also vector 
(shear) perturbations (without vorticity [14]), because it does not reduce to the form of Dijx in general. 
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III. CMB ANISOTROPIES DUE TO FIRST-ORDER AND SECOND-ORDER PERTURBATIONS 

In the unperturbed model universe, the observed temperature (To) of the CMB radiation is related to the 
emitted temperature (T e ) at the decoupling epoch (z e ) by T a = T e /(1 + z e ), and represented also as 

T = (uj /uj e )T e (3.1) 

using the observed and emitted frequencies w and w e (= (1 + z e )uj ). 

In the perturbed universe, these temperatures depend on the motions of matter and observers and on light 
paths passing through the inhomogeneous matter, and are expressed as 

T (x ,e) = (w /w e )T e (p,d) (3.2) 

with u) = —g^u^k", where is the observer's and emitter's velocities, k v {= dx u /dX) is the wave vector of 
photons with affine parameter A, (x Q , e) and (p, d) are (the position vectors and directional unit vectors) of the 
observer and emitter, respectively. 

The wave vector fc M satisfies the perturbed null geodesic equation and its solutions to the second-order are 
expressed as 

— k^ A- k>* A- k 1 * 

where x^ represents the light path and (r) denotes the r-order smallness. 
The temperature at the decoupling epoch is expressed as 

T e (p,d) =Tj°> [l+r(p,d)] (3.4) 

and the frequencies to the second-order are 



LO = 10 



(0)[i +(i (i) + £(2)]. (3.5) 



Then we have 



where 



T (x ,e) =T^[1 + 6T/T + ST/T], (3.6) 

1 2 



6 T/T = Cj^-uj^+t, 

5 T/T = u,™ - £jW + _ + _ CjW)t + pW&r/dx' + S^dr/dd\ (3.7) 

The procedure for solving null geodesic equations in perturbed universe models was shown by Pyne and Carroll, 
in which the background null rays are given by x^ ^ = (A, (A Q — X)e % ) and = (1, — e 4 ), and the baoundary 
conditions at the origin are x^^(X ) = x i - 2 ^{X ) = and k^ l (X ) = k^ l (X Q ) = 0. The expressions for Si T/T 
and 62 T/T were derived by Mollerach and Matarrese [7] in general gauges using Pyne and Carroll's procedure 
[13]. Their expressions are found to be gauge- invariant, and so the values can be calculated in a specified gauge 
without loss of any generality. Here we describe them using our solutions in the synchronous gauge, under the 

(r) (r) 

condition that the observers and emitters are comoving, i.e. v = Ve = 0. 

The following temperature perturbations are derived in the synchronous gauge from equations (2.20) - (2.28) 
in Mollerach and Matarrese's paper [7]. First we have the first-order perturbations: 

5T/T = r-h(X e ), (3.8) 

where 

h(X)=j dXA^ (A), (3.9) 

A (1) =^ (1) -^xg ) eV, (3.10) 
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and e l denotes a component of the directional unit vector e. The first-order wave vectors are 

fc«O(A) = -0« + ix^e i e j +7 1 (A) 



and 



where 



fc«*(A) = 20< 1 >e i - xptej - H (1) e l + X im ^ - J* (A), 
I{{\) = J dXA^ l {X) 



(3.11) 



(3.12) 



(3.13) 



and |i = 5 9; and ej — S^e 1 . 
The first-order light paths are 



z«°(A) = (A-A ) 

+ 



dX(X - X)A^'(X), 



P (i)« 



(A) = (A-AoJp^V-x^'e,-] 



- ^ rfA[2^ 1 )e i -x (1)ij e J - + (A-A)A( 1 )l i (A)]. 



The second-order wave vectors satisfy the following relation 

fc ( 2) o_ fc ( 2) o =/2(Ae)j 

where 



(3.14) 



(3.15) 



(3.16) 



I 2 (X) = J dX 



l -A^' + xf/e^k^ + e*fc«°) + 2k^A^' + 2<f>W A™ + x^ ^ 1 )" + z^^'l (3.17) 



and 



A^ = iv!?eV'. 



(3.18) 



The second-order temperature perturbations are 

= -/ 2 (A e ) + / 1 (A e )[/ 1 (A e )-r-^ 1 ) + i x ( 1 ^e 4 e i 
+ 4 1 )°4 1)/ + (4 1)j +4 1)0 e J -)r,, + g ? rfW, 

where 

d {1)l = e l - {e l - k {1 ^)/\e l - k {1)l \. 
Moreover, if we substitute our metric perturbations into the above equations, we obtain 



S T/T = t+- dXP'^Fijje* 



for the first-order perturbation. Since dP/dX = P' and dF/dX = —F ti e l , Eq.(3.21) can be expressed as 



5T/T = e 1 + e 2 



(3.19) 
(3.20) 

(3.21) 

(3.22) 



with 



©i = T--{(P'F, l ) e -(P'F^ ]e\ 
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e 



±J'd\P"(r 1 )F ti e i . (3.23) 



The latter term O2 represents the first-order Sachs- Wolfe effect. 
Next, we notice that Eqs.(3.11), (3.12), (3.14) and (3.15) lead to 



x (i)o + x (i)i e . = _ f dX AW, (3.24) 

JXo 



and we have 



f ' dXA^'(X)h(X) = hh(X e )} 2 . (3.25) 

J A G 

Then we obtain from Eq.(3.17) 

h{K) = \[h{K)] 2 (A e - Xo)AlVAW' + A?)' f° d\[AW + (A e - A)A«'] 

- 4 1 >/ 1 (A e ) + j\x{^A^' + A^A^' -AW" j\~XA^(X)}, (3.26) 
and therefore from Eq.(3.19) 

5T/T = h(X e )[^h(X e ) - t] -[A™' +r i e i ] f ' d\A™ 

^ Ae rfA {iA( 2 ) / + A( 1 )^ 1 )'-^ 1 ) / '^ A rfAA( 1 )(A)} + ^dW i , (3.27) 
where (ij, x l ) = (A, X D — X) in the integrands and 

/i(Ae) = ~\ J ' dXP'F l3 e l e\ 

A (2) = -^[PL{+P 2 Ml +QN ij + C 3 i ]e i e>. (3.28) 

In Eq.(3.27) the terms with path integrations represent the second-order nonlinear Integral Sachs- Wolfe effect, 
which brings the observational coupling of two linearly independent inhomogencitics with different wavelengths. 
As can be seen from Eq.(3.28), the induced gravitational radiation contributes to the CMB anisotropies. 

IV. PERTURBATIONS IN THE POISSON GAUGE 

In this section we derive the perturbations in the Poisson gauge which is defined by the condition 

z| r)l * = and = 0. (4.1) 

For this purpose, we use a gauge transformation from the perturbations in the synchronous gauge to those in 
Poisson gauge. The first-order gauge transformation has fully been studied by many authors (Bardeen [14], 
Kodama and Sasaki [15]). The second-order gauge transformation has more recently been derived by Bruni et 
al.[5] and the transformations of an arbitrary perturbed tensor T = % + ST^ + ^ST^ are expressed in terms 
of generators and £( 2 ) as 

5f{ 2) = 5T {2) +2 £ C(i)(5r (i) + £| (i) T + C im %, (4.2) 
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where C denotes the Lie derivative and the components of the generators are expressed as 

fO 

»(r) 



e M = « (r) (4.3) 



and 

C {r) = /? (r)|i + d {r)i (4.4) 

with d^\i = 0. 

This gauge transformation has been applied by Matarrese et al. to derive the second-order perturbations of 
the Einstein-de Sitter model in the Poisson gauge from those in the synchronous gauge. Here in the similar 
manner we will apply the transformation to our perturbations derived in §2. 

A. First-order transformation 

For the transformation of our metric perturbations from the synchronous gauge (S) to the Poisson gauge (P), 
we have 

a ' (i) 
a 











z (1) 




z Pi 






= W 








+ 2/3« -- 




= 0. 



Using Eq.(2.39) for (f>^ and Xs^"; we obtain from the above equations 

a« = --P'F, 
2 

/?« = ~PF, 



2 



St? IP = Spf/ P + P -a^\ 



P 

( p 1)0 = -^a«-a'«, 
a 



so that 



Similarly we use the transformations expressed as 

= Sg% + 2£ iw 8g% + g$ + C t(2) g<$ , 



(4.5) 



0. (4.6) 
The density and velocity perturbations satisfy the following relations: 



sp^/p = -L s (* P >-i\ AF+ i?Lp> Ft 

2pa z \ a J 2 a 

p 2 a\a J ' 

4 )l = \p'Fi- (4-8) 
B. Second-order transformation 
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5u (2)p = Su (2)s + ( c tm + £ |i ) ) u ( I o) + 2£ C(dKi)s- 



(4.9) 



More explicit expressions for these relations are given in Eqs. (5.8) - (5.14) of Matarrese et al.'s paper [6]. By 
analyzing these latter equations, we obtain the following generators 



where 



0* . p(fe„ + F=) + ,pT(fe + I^) + ^(-|*„ + lF,F I ) + | * . 



Ad 



(2) _ 



200 



0,i 



±(F >l F tl ) >i -F ti AF][p+±(P'y 



(4.10) 



200 



[F M F kl - (AF)% 



A9 = *o F t Fi. 



(4.11) 



It is found from Eq.(4.10) that vector perturbations without vorticity[14] appear also in this gauge. The resulting 
metric perturbations are expressed as 



50 r 6a' , 
i- — 2 P' + 



3 I a 



(4.12) 



(4.13) 



and 



(Azf } ) P = P'(l + P") [-^°*o,i + ^(WO.i - F,iAF 



(Xg }ll )p = and ^')p = Ci : 



(4.14) 



(4.15) 



where we used Eq.(2.8) in the derivation of the above equations, and || and T denote the scalar perturbation 
and the transverse and traceless part representing induced gravitational radiation, respectively. In the limit of 
A = 0, these generators and solutions reduce to Eqs. (6. 6), (6.7) and (6.8) in the Matarrese et al.'s paper [6], 
which were shown for the Einstein-de Sitter model. 
The density and velocity perturbations are 

^p/p) p = -ioo^Wi --^') e ° + 1 ?-(^'-^')^ + { 3 [(-) 2 - \~] ^ 

a V a ) 7 a \ 6/ I LV a / 4 oJ 



_ ~P'(P" + 2)}F 2 + ± [-2^PP' + (1 - ^P')/(pa 2 )]FjFj 

+ 4i - \P' \(- - 5(-) 2 )^' + 3-1 + 4(l - a -P') }FAF 
pa 1 I 2 LV a a / a J V a J i 

- \-PP'\3tAFf + 4F.uF.ul + rF|(AF) 2 + F„F >U ] 

( CI Z 

+ \(P'f[{AFf-F M F M ]}, 



(4.16) 
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/ n x 100 r a' /a'\2 / -i 100 r, /x9 / a' a 

K 8) )p = -_[i-3-j" + 2(-) (P') 2 ]e + ^[(P') 2 + P(i--P') 



+ 



6 V a 

(P') 2 --P(l--P' 
V a 



-l + ^P' + (^_ 2 (-) 2 )(P') i 
a V a a / 



(4.17) 



K 2) )p = -* {2) ; -^\1 + P")Q Q ,+ ™(PP' -\q')^o, 

- ^PP'FiFu-^P'i^-^P^FFi. (4.18) 

In the limit of A — 0, these relations also are consistent with Eqs. (6.10), (6.11) and (6.12) in Matarrese et al.'s 
paper [6], except for a few terms which may include some misprints with respect to (5^p)p and (u/ 2 ))f- 

Now we consider the second-order rotational velocity in the Poisson gauge. The rotational velocity vector 
(wp f ^) and the corresponding scalar quantity (wp') in the Poisson gauge are related to those (S) in the 
synchronous gauge by Eq.(4.2) as ST^ — Wp r ',Wp' and ST^ — ujg^ r \w^\ respectively, for r — 1, 2. Here we 
have T° = and 5T^ = for r = 1, 2, as was shown in §2. It is found therefore from these relations that to P 

(r) 

and ujp' vanish. That is, no rotational perturbations are induced also in the Poisson gauge. However vector 

(2) 

(shear) perturbations without vorticity appear in the form of Ad- ' . 



V. CONCLUDING REMARKS 



We have studied the second-order perturbations in spatially flat, pressureless cosmological models with 
nonzero cosmological constant, and could describe them explicitly using an arbitrary potential function F 
of spatial coordinates. It seems however to be difficult to derive similar results in the other cases with finite 
pressures. 

It is found in the second-order that tensor (gravitational- wave) perturbations and vector (shear) perturbations 
without vorticity are induced from the first-order scalar perturbations. 

Next we have derived the second-order temperature anisotropies of CMB in terms of these perturbations. 
They will be useful to analyze the nonlinear influence of local inhomogeneitics (at later stages) upon observed 
CMB anisotropies. Since nonlinearity brings the coupling of two linearly independent inhomogeneities with 
different wavelengths, it may appear as a small directional asymmetry and some deviations from the results 
which are expected in the standard linear theory. 
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APPENDIX A: SECOND-ORDER PERTURBED RICCI TENSORS AND EINSTEIN EQUATIONS 

We show here the expressions for perturbations of Ricci tensors to derive the perturbed Einstein equations. 
For the metric perturbations 5g^ v = h^ v + 1^ , the contravariant metric perturbations ig^ are derived from 
the condition {g^ v + Sg^ v )(g^ + 5g vl ) = 6^ as 

5g^ = -h» v + {h»W v -£» v ). (Al) 

The perturbed Christoflel symbols are 

= l^+^u-^-lKih^ + h^-K^ (A2) 
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where a semicolon denotes the four-dimensional covariant derivative in the unperturbed universe. From Eqs. 
(Al) and (A2), we can derive the perturbed curvature tensors and Ricci tensors: 

5R„ 7 — 8R^ IJ/ j = —(h^. rf ^ + h^ lv ^ — h 1/ j'. l ^ — h ;l/7 ), 

r DM If/"* -I-/"* — /> W — (>* — J- P J^M 

t^^vX-f — 2^ !y !7 A 7;' yA tly 7;A ^A;^ t i^;A7 ' l "*;7i 

— ^ai^-v-riX + ^";!/A — ^7;A _ ^A;j/7 — ^";A7 + ^Aj") 

+ l(fc£. 7 - + h ia ^){K. x + h% v - h^ a ) 

2 2^ 

2 a\ v:~{^ ' 7;^M v *1\\l "'p.-.u-yj 

1 



- i (2^ ;Al -/ i;Q )(/ i « 7 + / i « ;iy -/i I , 7 ;Q ). (A3) 



Using the relation 



+ + £i£ = (.g^ 7 + 8g^){R vl + 5 R v ~ f + 8 R„ 7 ) , (A4) 

we obtain the mixed components of the perturbed Ricci tensors: 

- g^SR v -y-h^R VJ , 
i i 

= SR.y - SR.j + (h^ ~ £^)R UJ . (A5) 

Here let us impose the synchronous gauge conditions (2.6) and (2.15) on /i M „ and Then we obtain the 

following expressions for second-order perturbed Ricci tensors S2 R%- 

On' n' 

2a 2 SRI = ^ + (il)" + ^(tiy + ySi-[hi(h^" 



2a 2 sR° = (£\iy-(^ k y+hM lk y-(h l kli y) 

l(hfyh l kli + 1 -(h^'(2h l kll -h lk ), 

2a 2 sRi = (&y -(ify + h kl (K u y -hUh^y -hhtyh 1 ? 



"k\V 'H V'fc ) 2 

+ \{K)'{2hf l - h \ k ) + K[{h\ k y -{hi !')'], 

2a 2 |i? ° = £" + ^'-^[(4)" + ^(4)']- ^(tf)W> 

2a 2 «J R = 2a 2 5(R°o + R k k ) = $1 + 2(£" + 3-£') 
22 a 

- [2 W + ) W - \{h'f + ^-hHhin (A6) 

where \i means the three-dimensional covariant derivative in the constant-curvature space with metric 7^, 
h = h k k , I = £ k and $^ is defined as follows: 



/ U-i 

k\i J 
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l^-k^+^h^-hf 1 ), (A7) 



where K (= ±1, 0) is the three-dimensional spatial curvature. In the Appendix A, we take the spatial curvature 
and the pressure p into consideration. 

The first-order and second-order perturbed Einstein equations are expressed as 

6K - \k{R = st^ 

SK - bt5R = 5T». (A8) 

2 2 2 2 

The perturbations of the energy-momentum tensor are 

ST* = 5ZSp + g„ a [u a u»(5p + 5p) + {5u a u» + u a 6u»)( P + p)} 

+ h ua u a u^{p + p), 
8T» = S^dp + g va [u a u^(5p + Sp) + (5uV + u a Su^)(Sp + Sp) 

2 2 221 111 

+ (6u a u^ + u a 6u^ + 6u a 6u^)(p + p)} 

+ h va [u a u»(S P + 5p) + (S uV + u a 5 u")(p + p)} 

+ t va u a ui>(p + p). 1 (A9) 



For velocity perturbations, we obtain 



using the identity 



Su° = 0, Su° = ^-g k i5u k 8u\ (A10) 
l 2 2a l l 



(.9^ + V + V)K + £ «" + *«")(«" + Su v + Su v ) = -1 (All) 



and the synchronous gauge conditions. The components of the energy-momentum tensor are, therefore, 



and 



ST?=5j6p, ST$ = -a{p + p)5u j , 5T° = -5p, (A12) 



ST? = 6j5p+(p + p)gi k 5u k 5u j , 

ST° = -g lk [Su k (Sp + Sp)]+Su k (p + p) + -h lk Su k (p + p), 

2 Cb 1 1 1 2 Qjl 

<5TJ = -aSu j (Sp + Sp) -adu j (p + p), 

2 111 2 

<5T ° = -Sp-(p + p)g k iSu k Su l . (A13) 

2 2 11 



If we assume the equation of state p = p(p), then we have 

1 



Sp= (dp/dp) 5 p and Sp=(dp/dp)Sp+^(d 2 p/dp 2 )(Sp) 2 , (A14) 



so that 

52? = -5i(dp/dp)5T°, 



i 

52? = -5/(dp/rfp)5T ° + i5f(rfV^ 2 )(^) 2 

2 2^1 



+ ip + p) [gik S u ] - 81 {dp/dp)g M S u l ] S u k . (A15) 
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When the first-order quantities h^SiU 1 , and Si p are known, the equations for determining the second-order 
metric perturbations are derived from Eqs.(A8) and (A15) as follows: 



2a' 

+ —h{(h^y + 2a 2 ( P + p)g tk 5u k S^, (t ^ j), 



(A16) 



+ 



In' In' 



1 



'(1 + dp/dp) + $• --$*(! + dp/dp) 
u. a z 

1 9o' 

htihfy - \^y{h{)' + \{ h 'f ^hUhfy 

2a 2 (p + p) g lk Su l - (dp/ dp)g k i Su l 5u k 
L 1 1 J 1 

{dp/d P )[\{h k )'{h l k y + ^hi(h k y (h'f] 

a 2 (d 2 p/dp 2 )(Sp) 2 , 



(A17) 



where in the latter equation the index i is non-dummy but others are dummy, h — h k , I = l\, and $^ is defined 
in Eq.(A7). Eqs.(A16) and (A17) consist of six equations, which corresponds to six independent components of 

4- 

The second-order density and velocity perturbations can be expressed in terms of metric perturbations as 



Sp/P 
2 



2pa 2 



5 Rl + \ 5 R - [p + p)gu 8u k 8u l 
2 22 11. 

2 ^t' -\(h k y(h[y - 2 ^h k (h[y + 1 -(hr 

,fe r„J 



+ -<&l-2a\p + p)g k i5u k 5u l 

L 11 



(A18) 



aSu = 

2 p + p 



SRn + aSu l (5p + Sp) 



2 {P h^{ {eliy -^ y+h ^ ly 

hM'y - \{h k yhf + \{h\)'(2ht -h* 

+ /4[(/il fe )' - (h kll )'] + 2a 3 6u i (5p + 6p)}. 



(A19) 
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